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Introduction 

The study of the fundamental theory of the strong interaction - Quantum Chromo Dynamics 
(QCD) - in extreme conditions of temperature, density and small parton momentum fraction 
(low-;c) has attracted an increasing experimental and theoretical interest during the last 20 
years. Indeed, QCD is not only a quantum field theory with an extremely rich dynamical 
content (asymptotic freedom, infrared slavery, (approximate) chiral symmetry, non trivial 
vacuum topology, strong CP violation problem, Ua{1) axial-vector anomaly, ...) but also 
the only sector of the Standard Model whose full collective behaviour - phase diagram, 
phase transitions, thermalization of fundamental fields - is accessible to scrutiny in the 
laboratory. The study of the many-body dynamics of high-density QCD covers a vast range 
of fundamental physics problems (Fig. [I]): 

• Deconfinement and chiral symmetry restoration: Lattice QCD calculations [IJ predict 
a new form of matter at energy densities (well) above £„,, ~ 1 GeV/fm-^ consisting of 
an extended volume of deconfined and bare-mass quarks and gluons: the Quark Gluon 
Plasma (QGP) |l2l. The scrutiny of this new state of matter - equation-of-state (EoS), 
order of the phase transition, transport properties, etc. - promises to shed light on 
basic aspects of the strong interaction such as the nature of confinement, the mechanism 
of mass generation (chiral symmetry breaking, structure of the QCD vacuum) and 
hadronization, which still evade a thorough theoretical description [|21 due to their highly 
non-perturbative nature. 
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Figure 1. Left: The multiple facets of QCD (adapted from H). Right: QCD phase 
diagram in the temperature vs. baryochemical potential (T.^b) plane. The arrows indicate 
the expected crossing through the deconfinement transition during the expansion phase in 
heavy-ion collisions at different accelerators. The dashed freeze-out curve indicates where 
hadrochemical equilibrium is attained in the latest stage of the collision ||5l. 



• Early universe cosmology: The quark-hadron phase transition took place some 10 ^us 
after the Big-Bang and is believed to have been the most important event in the Universe 
between the electro-weak (or SUSY) transition (x ~ lO^^'^ s) and nucleosynthesis (x ~ 
200 s). Depending on the order of the transitiorQ, several cosmological implications have 
been postulated [6] such as the formation of strangelets and cold dark-matter (WIMP) 
clumps, or baryon fluctuations leading to inhomogeneous nucleosynthesis. 

• Parton structure and evolution at small-x: HERA data flVl indicate that when probed 
at high energies, hadrons consist of a very dense system of gluons with small (Bjorken) 
momentum x = Ppa,,„„/ p,,adwn- At low x, the probability to emit an extra gluon is large, 
proportional to a^ln(l/;c), and gg fusion processes will play an increasing role in the 
parton evolution in the hadronic wavefunctions. At high virtualities and moderately 
low X, such evolution is described by linear DGLAP [[8l or BFKL [9J equations, suitable 
for a dilute parton regime. At ;c < 10^^ and below an energy-dependent "saturation 
momentum" Qs, hadrons are however more appropriately described as dense, saturated 
parton systems in the context of the "Colour Glass Condensate" (CGC) HIOI effective 
theory with the corresponding non-linear BK/JIMWLK [il2l [TTl evolution equations. 
Since the growth of the gluon density depends on the transverse size of the hadron, 
saturation effects are expected to set in earlier for ultrarelativistic heavy nuclei (for which 
2y oc A^/^, with A the number of nucleons) than for free nucleons. 

• Gauge/String duality: Theoretical applications of the Anti-de-Sitter/Conformal-Field- 
Theory (AdS/CFT) correspondence provide results in strongly coupled (i.e. large 't Hooft 
coupling X = g^ 1) SU(A'c) gauge theories in terms of a weakly-coupled dual gravity 

t The order itself is not exactly known: the transition, which is 1 *'-order in pure SU(3) gluodynamics and of a 
fast cross-over type for Nf = 2+1 quarks H], is still sensitive to lattice extrapolations to the continuum limit. 
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theory [fT3l . Recent applications of this fomialism for QCD-like = 4 super Yang- 
Mills) theories have led to the determination of transport properties of experimental 
relevance - such as the QGP viscosity [[T4|. the "jet quenching" parameter {q) [fT5l . or the 
heavy-quark diffusion coefficient [fT6l - from black hole thermodynamics calculations. 
Such results provide valuable insights on dynamical properties of strongly-coupled 
QCD that cannot be directly treated by perturbative or lattice methods, and open new 
phenomenological and experimental leads. 

• Compact object astrophysics: At high baryon densities and not too high temperatures, 
the attractive force between (colour antisymmetric) quarks can lead to the formation of 
bound {qq) Cooper pairs. Cold dense matter is thus expected to behave as a colour 
super-conductor with a non trivial quark pairing structure due to the combination of the 
various quantum numbers involved (spin, colour, flavour) [fTTl . This regime, currently 
beyond the direct reach of accelerator-based research (except indirectly in the region of 
baryon densities around the QCD tri-critical point. Fig. [T] right), may be realised in the 
core of compact (neutron, hybrid or other exotic) stars and, thus, open to study through 
astronomical observation. 

The only experimental means available so far to investigate the (thermo)dynamics of a multi- 
parton system involves the use of large atomic nuclei collided at ultrarelativistic energies. 
Figure [21 left, shows the total center-of-mass energy available for particle production (i.e. 
subtracting the rest mass of the colliding hadrons) at different accelerators as a function of 
the first operation year ("Livingston plot") [18]. The exponential increase in performance 
translates into an energy doubling every 2 (3) years for the ion (p,p) beams. Head- 
on collisions of heavy ions (AA) can produce extremely hot and dense QCD matter by 
concentrating a substantial amount of energy (<^(1 TeV) at mid-rapidities at the LHC, see 
Fig. |2] right) in an extended cylindrical volume V — tiR^Xq ~ 150 fm^ for a typical large 
nucleus with radius Ra = 6.5 fm, at thermalization times of Xq = 1 fm/c. 




Figure 2. Left: "Livingston plot" for (anti)proton and ion accelerators in the period 1960- 
2008 (adapted from ifTSl ). Right: Measured transverse energy per unit rapidity at r| = 0, and 
corresponding Bjorken energy density eB,(To = 1 fm/c) |fT9l . in central heavy-ion collisions at 
various cm. energies Ii20i i2 1 j fitted to a logarithmic parametrization. 
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The hot and dense systems produced in high-energy AA collisions are not prepared under con- 
trolled thermodynamical conditions but they follow a dynamical trajectory along the phase 
diagram shown e.g. in Fig. [T} right. After the collision, the system (with a temperature 
profile decreasing from the center) expands with relativistic longitudinal (transverse) veloc- 
ities (P) ~ 1.0(0.5) and cools at rates T oc x^^l" (e.g. n = 3 for a longitudinal-only expan- 
sion [19]). When T reaches 7^.,,, ~ 190 MeV, the quark matter undergoes a phase transi- 
tion into hadrons. The produced hadronic gas stops self-interacting collectively at freeze-out 
times X ~ 10-20 fm/c f2^. At the initial stages of the reaction (1 fm/c after impact), the 
commonly used "Bjorken estimate" [19| gives energy densities attained at mid-rapidity of 
Efi, = dET/dV[\^=Q/{nR\xQ) ^ 5, 10 GeV/fm^ at RHIC and LHC (Fig. E] right). Although 
these values can only be considered as a lower limit since they are obtained in a simple 1-i-lD 
expansion scenario ignoring any effects from longitudinal work, they are already about 5 and 
10 times larger, respectively, than the QCD critical energy density for deconfinement. High- 
energy collisions in heavy-ion colliders provide therefore the appropriate conditions for the 
study of highly excited quark-gluon matter. 

This review is organised as follows. Sections \T\ introduces the experimental probes of 
subhadronic matter used in high-energy AA collisions. Section [2] briefly reviews the 
experimental apparatus used in heavy-ion collisions at RHIC and SPS. In Section [3] various 
hard QCD results from proton-proton collisions at RHIC are presented as the "free space" 
baseline to which one compares the heavy-ion ("QCD medium") data. Sections \4\WT\ each 
discuss a different physics observable in the context of the latest experimental results available 
at RHIC and SPS. Due to space limitations, the list of topics chosen covers but a fraction of the 
substantial amount of data collected in the last years. Those observables that provide direct 
information on the partonic phases of the reaction have been given preference over "soft" or 
bulk observables from the late hadronic stages. Thus, important results relative to the freeze- 
out phase of the collision - e.g. chemical and kinetic equilibrium from light (and especially 
strange [|23ll24l ) hadron abundances [25], "femtoscopy" measurements from HBT radii [|26ll . 
or possible signatures of (prehadronic) divergent susceptibilities from final net charge and 
{Pt) fluctuations [|27l - are not treated here. The selected experimental results are mostly 
from the comprehensive reviews of the four RHIC experiments (PHENIX [12811 . STAR [i29ll . 
PHOBOS [30J, and BRAHMS [i31J ) from AmAm, dAu and pp collisions up to a maximum 
center-of-mass energy of = 200 GeV, as well as recent results from NA60 Inin reactions 
at SPS ( ^ = 17.3 GeV) [|32l. 

1. Experimental probes of QCD matter 

Direct information on the thermodynamical and transport properties of the strongly interacting 
medium produced in AA collisions is commonly obtained by comparing the results for a 
given observable ^aa to those measured in proton(deuteron)-nucleus {p{d)A, "cold QCD 
matter") and in proton-proton (pp, "QCD vacuum") collisions as a function of cm. energy, 
transverse momentum, rapidity, reaction centrality (impact parameter b), and particle type 
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The observed enhancements and/or suppressions in the RAA,dA{ \/s^,PT,y,tn',b) ratios can 
then be directly linked to the properties of the strongly interacting matter after accounting for 
a realistic modeling of the space-time evolution of the AA expansion process. 



Among the observables whose suppression is expected to provide information on the produced 
system, we will discuss: 

• The total particle multiplicity which, related via local parton-hadron duality [|33l to 
the number of initially produced partons, will be suppressed if the initial parton flux is 
reduced due to low-x saturation effects in the colliding nuclei L34J (Section H]). 

• High-pr leading hadrons are expected to be produced in reduced yields due to medium- 
induced energy loss via gluonsstrahlung of the parent partons in a system with a large 
number density of colour charges JSSllMl (Sections |6] and |7]). 

• Dissociation of the heavy quarkonia bound states has long since proposed [|37l as a 
sensitive signature of Debye screening effects above 7^,„, (Section [TOl). 

• The mass of light vector mesons is expected to drop in some scenarios of chiral 
symmetry restoration ll38l [39l which directly link the in-medium meson mass to a 
temperature- (and baryon density-) dependent chiral condensate: my my ■ {qq) [T) 
(Section [n])- 

Likewise, the following observables discussed hereafter are expected to be enhanced in a 
strongly-interacting multi-parton system compared to the results m pp collisions: 

• The soft hadron spectra {pT < 2 GeV/c), both inclusive and relative to the reaction 
plane, are expected to be boosted due to the development of collective radial and elliptical 
(hydrodynamical) flows in the early partonic phases of the reaction [|22ll40l (Section [5]). 

• Semihard {pT ^ 2-5 GeV/c) yields and flows of baryons are predicted to be enhanced in 
the context of parton recombination (or coalescence) models [ 41| which take into account 
modifications of the hadronization process in a dense deconfined medium (Section[8]). 

• The measurement in AA collisions of a thermal photon excess over the expected prompt 
Y yield ||42|| would provide direct access to the thermodynamical properties (temperature, 
EoS) of the produced system (Section HI). 

Among all available experimental observables, particles with large transverse momentum 
Pt and/or high mass ("hard probes") [|43l l44l are of crucial importance for several reasons 
(Fig. [3l): (i) they originate from partonic scatterings with large momentum transfer 
and thus are directly coupled to the fundamental QCD degrees of freedom; (ii) their 
production time-scale is very short, x I / pr '^O-l fm/c, allowing them to propagate through 
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Figure 3. Examples of hard probes whose modifications in high-energy colHsions provide 
direct information on properties of QCD matter such as the {q) transport coefficient, the initial 
gluon rapidity density dN^ /dy, and the critical temperature and energy density. 



(and be potentially affected by) the medium, (iii) their cross sections can be theoretically 
predicted using the perturbative QCD (pQCD) framework. Hard processes thus constitute 
experimentally- and theoretically-controlled self-generated "tomographic" probes of the 
hottest and densest phases of the reaction. The pQCD factorization theorem B31 allows one 
to determine the production cross section of a given hard probe as the convolution of long- 
distance parton distribution (PDFs, fa/A) fragmentation (FFs, Dc^h) functions and the 
(perturbatively computable up to a given order in a^) parton-parton scattering cross section: 

d<^ABih = fa/A{x, Q^)®fb/B{x, Q^) ® do'^^tc® D,^h{z, Q") + ^(l/Q^) (3) 

The validity of Eq. (|3]) holds on the possibility to separate long- and short-distance effects with 
independent QCD time (length) scales as well as on the assumption of incoherent parton- 
parton scatterings. In AA collisions, the incoherence condition for hard processes implies: 
fajA ~ ^ ■ /a/A'' i-^- that, in the absence of medium effects, the parton flux in a nucleus A 
should be the same as that of a superposition of A independent nucleons. Thus, 

do%th-^-'^ ■ falpM^)® fblp{x,Q')®d<5'2tc®D,^h{z.Q'), (4) 

and minimum-bias hard cross sections in AB collisions are expected to scale simply as 
do^'^\^, =A B- do'^f. In the most general case, for a given AB reaction with arbitrary 
impact parameter b the yield can be obtained by multiplying the cross sections measured 
m pp collisions with the ratio of the incident parton flux of the two nuclei: dN^^g'^{b) — 
{TAsib)) ■d(5^p'^, where TAsip) (normalised to A • B) is the nuclear overlap function at b 
determined within a purely geometric Glauber eikonal model using the measured Woods- 
Saxon distribution for the colliding nuclei P6l . The standard method to quantify the effects 
of the medium on the yield of a hard probe produced in a AA reaction is thus given by the 
nuclear modification factor. 

V ( u\ d^NAA/dydpT 

RAA{pT,y;b) = ^ — , (5) 

{TAA{b)) X d^Opp dydpT 
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which measures the deviation of AA at b from an incoherent superposition of A^A'^ collisions. 
This normalization is usually known as "binary collision scaling". 

2. Experiments in high-energy heavy-ions physics 

The Relativistic Heavy Ion Collider (RHIC) pTj at Brookhaven National Laboratory is a 
3.8-km circumference accelerator composed of two identical, quasi-circular rings of super- 
conducting magnets (~400 dipoles and ~500 quadrupoles) with six crossing-points. The 
machine, which started operation in 1999, can accelerate nuclei (protons) up to a maximum 
of 100 (250) GeV/c per nucleon. The center-of-mass energies in AA and p{d)A collisions, 
^/s^ = 200 GeV ( = 500 GeV for pp), are more than an order of magnitude larger than 
those at the CERN SPS i^/s^= 17.3 GeV). The currently attained average A4 luminositie^ 
(=Sf) ^ 4 X 10^^ cm^^s^^ = 0.4 mb^^s^^ are twice the design luminosity thanks mainly to 
an improved vacuum system [48J. There are four dedicated experiments at RHIC: two large 
multi-detector systems (PHENIX and STAR, Fig. H]) and two smaller specialised spectrome- 
tres (BRAHMS and PHOBOS). 

BRAHMS P9l has two movable magnetic spectrometer arms with hadron identification 
(n,K,p) capabilities up to very large rapidities (ymax ~ 4 for charged pions). The two spec- 
trometers consist of a total of 5 tracking chambers and 5 dipole magnets (maximum field of 
1.7 T), plus 2 Time-Of-Flight systems and a Ring Imaging Cerenkov Detector (RICH) for 
particle identification. The angular coverage of the forward spectrometer (FS) goes from 2.3° 
< < 30° (solid angle of 0.8 msr) up to momenta of 35 GeV/c. The midrapidity spectrometer 
(MSR) covers 30° < < 95° (acceptance 6.5 msr). The PHOBOS experiment f50\ is based 
on silicon pad detectors and covers nearly the full solid angle (11 units of pseudorapidity) for 
charged particles, featuring excellent global event characterization at RHIC. PHOBOS con- 
sists of 4 subsystems: a multiplicity array ("octagon", |ri| <3.2, and "rings", |ri| <5.4), a finely 
segmented vertex detector, a two-arm magnetic spectrometer (with dipole magnet of strength 
1.5 Tm) at midrapidity including a time-of-flight wall, and several trigger detectors. Particle 
identification (PID) is based on time-of-flight and energy loss in the silicon. 

PHENIX [|5TI is a high interaction-rate experiment specifically designed to measure hard 
QCD probes such as high-pr hadrons, direct y, lepton pairs, and heavy flavour. PHENIX 
(Fig. m left) achieves good mass and PID resolution, and small granularity by combining 
13 detector subsystems (~350,000 channels) divided into: (i) 2 central arm spectrometers for 
electron, photon and hadron measurement at mid-rapidity (|ti | < 0.35, A(|) = tt); (ii) 2 forward- 
backward spectrometers for muon detection (|ri| = 1.15 - 2.25, A(|) = 2%); and (iii) 4 global 
(inner) detectors for trigger and centrality selection. Two types of electromagnetic calorime- 
ters, PbSc and PbGl, measure y and e^. Additional electron identification is possible thanks to 
the RICH detector. Charged hadrons are measured in the axial central magnetic field (strength 

§ Note that the "equivalent-/5/5 luminosity" for hard processes, obtained scaUng by the number of nucleon- 
nucleon collisions, is a much larger value: {^) pp^equiv ■ {■^)aa ~ /jb^'s^'. 



Quark-Gluon Matter 



8 




Figure 4. The two large experiments at RHIC: PHENIX IST] (left) and STAR 1521 (right). 



1.15 Tm) by a drift chamber (DC) and 3 layers of MWPC's with pad readout (PC). Hadron 
identification (tt^, K^, and p,p) is achieved by matching the reconstructed tracks to hits in 
a time-of-flight wall (TOF). Triggering is based on two Beam-Beam Counters (BBC, |ri| = 
3.0 - 3.9) and the Zero Degree Calorimeters (ZDC, with |0| < 2 mrad). PHENIX features 
a state-of-the-art DAQ system capable of recording ~300 MB/s to disk with event sizes of 
~100 KB (i.e. coping with event rates of ~3 kHz). 

The STAR experiment [[52l is based around a large-acceptance Time Projection Chamber 
(TPC) inside a solenoidal magnet with radius 260 cm and maximum field strength 0.5 T 
(Fig. HI right). The TPC with radius 200 cm and full azimuthal acceptance over |ri| <1.4 pro- 
vides exceptional charged particle tracking and PID via ionization energy loss dE jdx in the 
TPC gas and reconstruction of secondary vertices for weakly decaying particles. Additional 
tracking is provided by inner silicon drift detectors at midrapidity and forward radial-drift 
TPCs at 2.5< |ri| <4. Photons and electrons are measured in the Barrel and Endcap Electro- 
magnetic Calorimeters (EMC) covering -1.0< T] <2.0 and full (|). The large STAR coverage 
allows for multi-particle correlation studies, jet reconstruction in pp, and also measurement of 
strange and charm hadrons. Triggering is done with the ZDCs, forward scintillators (Beam- 
Beam counters), a barrel of scintillator slats surrounding the TPC, and the EMC. 

NA60 ll53l is a high interaction-rate ^xeJ-tor^e? experiment focused on the study of dimuon, 
vector meson and open charm production in pA and AA collisions at the CERN SPS. The 
17 m long muon spectrometer previously used by NA38 and NA50, is composed of 8 multi- 
wire proportional chambers, 4 scintillator trigger hodoscopes, and a toroidal magnet. This 
spectrometer is separated by a ~ 5 m long hadron absorber (mostly carbon) from the vertex 
region, where a silicon tracker made of pixel detectors ll54l placed in a 2.5 T dipole field, 
measures the produced charged particles. The matching between the muons and the vertex 
tracks leads to an excellent dimuon mass resolution. 



Quark-Gluon Matter 



9 



3. Benchmark channels in the vacuum and in cold QCD matter 

Proton-proton collisions are the baseline free-space reference to which one compares the 
AA results in order to identify initial- or final- state effects which modify the expectation 
of equation dH) and can thus provide direct information on the underlying QCD medium. 
Figure [5] collects six different /^j- -differential inclusive cross sections measured at RHIC in 
pp collisions at ^/s = 200 GeV: jets [|55l. charged hadrons ll56l . neutral pions [[57]|. direct 
photons [58], D,B mesons (indirectly measured via from their semileptonic decay) ||59l 
all at central rapidities (y = 0), and negative hadrons at forward rapidities (Tj = 3.2) [|60l . The 
measurements cover 9 orders of magnitude in cross section (from 10 mb down to 10 pb), 
and broad ranges in transverse momentum (from zero momentum for D, B mesons up to ~45 
GeV/c, a half of the kinematical limit, for jets) and rapidity (from r| = up to T] = 3.2 for h~ 
and even, not shown here, y = 3.8 for 7l°'s tMJ)- Whenever there is a concurrent measurement 
of the same observable by two or more RHIC experiments, the data are consistent with each 
other (the only exception being the heavy-flavour electron cross sections, which are a factor 
2-3 larger in STAR il than in PHENIX ll59l ). 




Figure 5. Left: Compilation of hard cross sections in pp collisions at = 200 GeV 
(10%-30% syst. uncertainties not shown for clarity purposes) measured by STAR 1551 156]| . 
PHENIX f57l|58l|59l, and BRAHMS |60| compared to NLO |63, 64| and NLL |65l pQCD 
predictions (yellow bands). Right: Nuclear modification factors for J/\\t versus rapidity ItTI 
(top), and for high-z^T- K° at mid-rapidity ll72l (bottom) in central dAu collisions at ^^i^J^ = 200 
GeV compared to pQCD calculations iClllTll with EKS98 ITS) nuclear shadowing. 

Standard next-to-leading-order (NLO) [t 63l [64ll or resummed next-to-leading log (NLL) |[65l 
perturbative QCD calculations with modern proton PDFs [66l, hadron fragmentation func- 
tions [[67l[68l . and with varying factorization-renormalization scales {/j = pt/2 — Ipj) show 
an overall good agreement with the available pp data at a/s = 200 GeV (yellow bands in Fig. 
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[51 left). This is true even in the semi-hard range />r ~ 1 — 4 GeV/c, where a perturbative de- 
scription would be expected to give a poorer description of the spectra. These results indicate 
that the hard QCD cross sections at RHIC energies are well under control both experimen- 
tally and theoretically in their full kinematic domain. This is at variance with measurements 
at lower (fixed-target) energies where several data-theory discrepancies still remain [169117011 . 

Not only the proton-proton hard cross sections are well under control theoretically at RHIC 
but the hard yields in deuteron-nucleus collisions do not show any significant deviation 
from the perturbative expectations. Figure [H right, shows the nuclear modification factors 
measured in dAu collisions at ^/s^ = 200 GeV for 7/\|/, RdAu{y) [EH, and for leading 
RdAuiPr) at J = [1721 . At mid-rapidity, the maximum deviation from the RdAu = 1 
expectation for hard processes without initial-state effects is of the order of ~20%. Both 
RdAuiy.Pr) ratios are well accounted for by standard pQCD calculations [[7311741 that include 
DGLAP-based parametrizations of nuclear shadowing i751 and/or a mild amount of initial- 
state pt broadening [i76ll to account for a modest "Cronin enhancement" \J7j. These data 
clearly confirm that at mid rapidities, the parton flux of the incident gold nucleus can be 
basically obtained by geometric superposition of the nucleon PDFs, and that the nuclear 
{x, Q^) modifications of the PDFs are very modest. Since no final-state dense and hot system 
is expected to be produced in dAu collisions, such results indicate that any deviations from 
Raa = 1 larger than ~40% potentially observed for hard probes in AuAu collisions (at central 
rapidities) can only be due to final-state effects in the medium produced in the latter reactions. 

4. Low-.x: gluon saturation: AA rapidity densities, and high-p^ dA forward suppression 

The bulk hadron multiplicities measured at mid-rapidity in central AuAu at ^/s^ = 200 GeV 
dNch/dr[\^=o ^ 700, are comparatively lower than the dNch/dr{\r\=o ^ 1000 expectations 
of "minijet" dominated scenarios [TTSl . soft Regge models [[79l (without accounting for 
strong shadowing effects [[80l ). or extrapolations from an incoherent sum of proton-proton 
collisions MM (Fig.[6l left). However, Colour Glass Condensate (CGC) approaches [[34l[82l 
which effectively take into account a reduced initial number of scattering centers in the nuclear 
PDFs, fa/Ai^iQ^) < ^ ■ /a/A'('^5 2^)' agree well with experimental data. In the saturation 
models non-linear effects become important and start to saturate the parton densities when the 
area occupied by the partons becomes similar to that of the hadron, nR^. For a nucleus with 
A nucleons and total gluon distribution xG{x, Q^) this condition translates into the following 
saturation momentum [[831 [84l : 

Q^(x) - a,-^xG{x,Q^) ~a1/3x-^ ~a1/3( V^)^ ~ A^/^e^^, (6) 

with X ~ 0.2-0.3 [[341 . The mass number dependence implies that, at comparable energies, 
non-linear effects will be A^^^ ^ 6 times larger in a heavy nucleus (A ~ 200 for Au or Pb) 
than in a proton. Based on the general expression ©, CGC -based models can describe the 
centrality and cm. energy dependences of the bulk A4 hadron production (Fig.[6l right). 
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Figure 6. Left: Data vs. models for dNch/dT\ \ty=Q in central AuAu at ./s^ = 200 GeV ll30l[8T1l 
(the saturation model prediction is identified as McLV |85|). Right: Energy- and centrality- 
(in terms of the number of nucleons participating in the collision, A^part) dependences of 
(iA^c/,/iiri|r|=o (normalised by A^part): PHOBOS data [301 versus the saturation prediction lf82l . 



The second manifestation of CGC-like effects in the RHIC data is the BRAHMS observation 
of suppressed yields of moderately high-pr hadrons {pT ~ 2 — 4 GeV/c) in dAu relative to pp 
collisions at forward rapidities (Tj 3.2, Fig.|7]) [|60l . Hadron production at such small angles 
is theoretically sensitive to partons in the Au nucleus with ^2'" = {pr / ^Aw) exp(— Tj) ^ 
^(10^^) fJM . The observed nuclear modification factor, RdAu ~ 0-8, cannot be reproduced 
by pQCD calculations that include the same leading-twist nuclear shadowing 11741 [86l [871 that 
describes the dAu data at 3^ = (Fig. [51 right) but can be described by CGC approaches that 
parametrise the Am nucleus as a saturated gluon wavefunction [[8811891 [90l . 




Figure 7. Nuclear modification factor RdAuipr) for negative hadrons measured at forward 
rapidities by BRAHMS in dAu at ^JJ^ = 200 GeV li60l compared to predictions of leading- 
twist shadowing pQCD |[74l[86l (two upper curves) and CGC |[88l[89l (lower curves). 

It is worth noting, however, that at RHIC energies the saturation momentum is in the transition 
between the soft and hard regimes {Q^ ~ 2 GeV^) and that the results consistent with the CGC 
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predictions are in a kinematic range with relatively low momentum scales: {pj) ~ 0.5 GeV 
for the bulk hadron multiplicities and {pj) ~ 2 GeV for forward inclusive hadron production. 
In such a kinematic range non-perturbative effects can blur a simple interpretation based on 
partonic degrees of freedom alone. Following Eq. Q, the relevance of low -x QCD effects will 
certainly be enhanced at the LHC due to the increased: center-of-mass energy, nuclear radius 
(A^/^), and rapidity of the produced partons f.76i,i91| . At the LHC, the saturation momentum 
2^ ^ 5 - 10 GeV^ |l34l will be more clearly perturbative and the relevant x values in AA and 
pA collisions will be 30-70 times lower than at RHIC: X2 ~ 10^^(10^^) at central (forward) 
rapidities for processes with a hard scale pr ~10 GeV. 

5. sQGP viscosity: Strong hydrodynamical flows and AdS/CFT connection 

The bulk hadron production {pj < 2 GeV/c) in AuAu reactions at RHIC shows strong 
collective effects known as radial and elliptic flows. First, the measured single hadron 
Pt spectra have an inverse slope parameter T^^ larger than that measured in pp collisions, 
increasing with reaction centrality and with hadron mass as expected if collective expansion 
blue- shifts the hadron spectra. Empirically, one finds: Zff^T + {^t) ' with T and (pr) 
being the freeze-out temperature and average collective flow velocity of the "fireball", and m 
the mass of the hadron. Phenomenological fits of the spectra to "blast wave" models yield 
transverse flow velocities (Pr) ~ 0.6 [28J. Full hydrodynamical calculations which start with 
a partonic phase very shortly after impact (Xo < 1 fm/c) develop the amount of collective 
radial flow needed to accurately reproduce all the measured hadron spectra (Fig. [8]). 
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Figure 8. Transverse momentum spectra for pions, kaons, and (anti)protons measured at 
RHIC below pj ~ 3 GeV/c in 0-10% central AmAm collisions at ^/s^ = 200 GeV compared 
to hydrodynamics(+pQCD) calculations [|92l . 
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Figure 9. Left: Spatial asymmetry with respect to the reaction plane of the produced "fireball" 
in non-central nucleus-nucleus collisions. Right: Measured elliptic flow parameter viipr) 
below pt = 2 GeV/c for a variety of hadrons 197] compared to hydrodynamic predictions ll40ll . 



Secondly, the azimuthal distributions dN / dA(^ of hadrons emitted relative to the reaction plane 
(A(|) = (|) — 4>^/)) show a strong harmonic modulation with a preferential "in-plane" emission in 
non-central collisions. Such an azimuthal flow pattern is a truly collective effect (absent in pp 
collisions) consistent with an efficient translation of the initial coordinate-space anisotropy 
in AA reactions with non-zero impact parameter (i.e. with a lens-shaped overlap zone, see 
Fig.|9]left) into a final "elliptical" asymmetry in momentum-space. Rescattering between the 
produced particles drives collective motion along the pressure gradient, which is larger for 
directions parallel to the smallest dimension of the lens. The strength of this asymmetry 
is quantified via the second Fourier coefficient V2(/'r,y) = {cos {2A(^)) of the azimuthal 
decomposition of single inclusive hadron spectra with respect to the reaction plane |l93l l94l 



The large V2 ~ 0.16 measured in the data (Fig. HI right) indicates a strong degree of collectiv- 
ity (pressure gradients) building up in the first instants of the collision. Indeed, elliptic flow 
develops in the initial phase of the reaction and quickly self-quenches beyond x ~ 5 fm/c as 
the original spatial eccentricity disappears El. Two additional experimental observations 
support the existence of an efficient hydrodynamical response with very short thermaliza- 
tion times. First, not only light hadrons but also D,B mesons (indirectly measured via their 
semileptonic decays into e^) show momentum anisotropics with V2 as large as 10% ||96l . The 
fact that the heavy c and b quarks participate in the common flow of the medium is clearly 
suggestive of a robust collective response during the early partonic phase. In addition, the 
V2 values measured for different centralities, at different center-of-mass energies (200 and 62 
GeV) and for different colliding systems {AuAu and CuCu) are found to show simple scaling 
laws with the reaction eccentricity |l98l HH also in agreement with hydrodynamics expecta- 
tions. 
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Figure 10. Elliptic flow (normalised by the participant eccentricity f98l) V2/8 as a function of 
the hadron rapidity density dNch/dr[ |r|=o normalised by the reaction overlap area Ax, measured 
at SPS [104] and RHIC ll98l 11051 compared to the "hydrodynamical hmit" expectations for a 
fully fliermaUsed system with hard (HRG-like) or soft (QGP-like) EoS [|95]|T0l[l08l. Adapted 
from Ill05ll98l (10% errors have been added to approx. account for V2, e syst. uncertainties). 



These results have attracted much attention for several reasons. First, the strong V2 seen 
in the data is not consistent with the much lower values, V2 < 6%, expected by transport 
models of hadronic matter HI 0011 or for a partonic system rescattering with perturbative cross 
sections (Ogg ~ 3 mb) HIOIL The magnitude, and the pr and hadron mass dependences 
of the radial and elliptic flows below pr ^ 2 GeV/c are, however, well described by ideal 
hydrodynamics models whose space-time evolution starts with a realistic QGP Equation of 
State (EoS) with initial energy densities £0 ~ 30 GeV/fm^ at thermalization times Xq ~ 0.6 
fm/c [I22I l40l 11021 1103II (Fig. |9l). Second, such a degree of accord between relativistic hy- 
drodynamics and the data was absent at lower CERN-SPS energies [il04il . Fig. [TO] shows 
the particle-density dependence of the V2 parameter (scaled by the eccentricity of the reac- 
tion £ to remove centrality-dependent geometrical effects) measured in semicentral collisions 
at different cm. energies. The fact that SPS data lie below the "hydrodynamical limit" 
curves [11081 l95l estimated for a system completely thermalised, suggests that equilibration 
is only partially achieved at the top SPS energy H 10411 . RHIC V2 data in the range ^/s^ ^ 62 
- 200 GeV [1105111061 l98l are, however, close to the full thermalization expectations. Third, 
inclusion of viscous (i.e. "internal dissipation") corrections to the ideal fluid dynamics equa- 
tions spoils the reproduction of V2{pt), especially above pr ^ ^ GeV/c where even a modest 
viscosity brings V2 towards zero [|109|[ . Estimates of the maximum amount of viscosity al- 
lowed by the V2{pt) data pi 101 give a value for the dimensionless viscosity/entropy ratio 
close to the conjectured universal lower bound, T\/s = fi/{An), obtained from AdS/CFT cal- 
culations [fT4| . Similarly, approaches [[111[[112[[ that describe simultaneously the large V2{pt) 
and quenching factors Raa{pt) of heavy-flavour e^ require small heavy-quark diffusion coef- 
ficients (3 < InT D < 6) which correspond to very small shear viscosities, 1 < 4n{r\/s) < 2, 
and/or very short thermalization times. 
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The fast (local) thermalization times supported by the robust collective flow generated in the 
first instants of the reaction, and the good agreement of the pt- and mass-differential spectra 
and V2 with ideal relativistic hydrodynamics models which assume a fluid evolution with zero 
viscosity (i.e. with negligible internal shear stress), have been presented as evidences that 
the matter formed at RHIC is a strongly interacting QGP (sQGP) hi f3l [TT4l [TT5l [TT6l fTTTll . 
This new state of matter - with liquid-like properties e.g. a Coulomb coupling parameter 
r = (Epot)/ {Eki„) ~ ^2(4i/3r)/(3r) ~ 3 for ^2 _ 4 - 6 at r ^ 200 MeV MEW - chal- 
lenges the anticipated paradigm [2J of a weakly interacting gas of relativistic partons (with 
r ^ 1), lending support to the application of strongly-coupled-gauge/weakly-coupled-gravity 
duality techniques [[TH [151 [T6]| to compute relevant sQGP parameters. 

It is worth noting, however, that recent lattice results predict a QCD transition temperature, 
Tcri, ~ 190 MeV [|119|1 , which is ~30-40 MeV larger than the freeze-out temperature extracted 
from observed particle yields in heavy ion experiments (dotted-dashed curve in Fig.[T]) Il5ll25l. 
Thus, an intermediate regime between the QCD transition and freeze-out could exist during 
which the system created in heavy ion collisions persists in a very dense strongly-interacting 
hadronic phase. At LHC energies the contribution from the QGP phase to the collective 
particle flow(s) will be much larger than at RHIC or SPS and, therefore, the V2 will be 
less dependent on the details of the subsequent hadronic phase. The measurement of the 
differential elliptic flow properties in AA collisions at the LHC will be of primary importance 
to confirm or not the sQGP interpretation as well as to search for a possible weakening of the 
V2 indicative of the existence of a weakly interacting QGP phase at higher temperatures than 
those of the Uquid-like state found at RHIC [QTSHTTOl . 

6. Parton number density and (q) transport coefficient: High-p^ hadron suppression 

Among the most exciting results of the RHIC physics programme is the observed strong 
suppression of high-pr leading hadron spectra in central AA [11201 consistent with the pre- 
dicted attenuation of the parent quark and gluon jets in a dense QCD medium ("jet quench- 
ing") [|121lll22n . Above pr ^ 5 GeV/c, neutral mesons (71°, Tj) [1123] and inclusive charged 
hadrons [[56l all show a common factor of ~5 suppression compared to an incoherent su- 
perposition of pp collisions (Fig. [TTl left). Such a significant suppression was not observed 
at SPS where, after reevaluating the pp baseline spectrum [|124[ \125L the central-AA me- 
son spectra show a R^a around unity (Fig. [HI right) probably due to the cancellation of a 
~50% final-state suppression by initial-state Cronin broadening [|124[ [1261 [TITl . At RHIC, 
the Raa ~ 1 perturbative expectation which holds for other hard probes such as "colour blind" 
photons ni28l - and for high-pj hadrons in dAu reactions (Fig. [51 right) - is badly broken 
{Raa ~ 0.2) in central AuAu collisions. This strongly supports the picture of partonic energy 
loss in final-state interactions within the dense matter produced in the reaction. 

The dominant contribution to the energy loss is believed to be of non-Abelian radiative nature 
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Figure 11. Nucleai- modification factor Raa{pt) for y, 7t°, r] ifTH . f (961, and /z± 1561 (left) 
and for "non-photonic" from D,B mesons |,62J (right) in central AuAu at y/s^ = 200 GeV 
compared to various parton energy loss model predictions II 1301 11331 11361 fTT2l . 



(i.e. due to gluon radiation) as described in the GLV lfT29l fBOll and BDMPS [|35l [BTTI 
(or LPCI II132II ') formalisms. In the GLV approach, the initial gluon density dN^ /dy of the 
expanding plasma (with transverse area and length L) can be estimated from the measured 
energy loss Afi: 

AEocalCR^^L. (8) 
dy 

where Cr is the Casimir colour factor of the parton (4/3 for quarks, 3 for gluons). In 
the BDMPS framework, the transport coefficienjjj] (q), characterizing the squared average 
momentum transfer of the hard parton per unit distance, can be derived from the average 
energy loss according to: 

{AE)ocasCR{q)L\ (9) 

From the general Eqs. ([8]) and very large initial gluon rapidity densities, dN^ /dy 1100 
± 300 IBOil, or equivalently, transport coefficients (<?) ^ 11 ± 3 GeV^/fm ifTM [T35l [T36l 
I137L are required in order to explain the observed amount of hadron suppression at RHIC. 
The corresponding values for SPS are dm /dy ^ 400 ± 100 and {q)K3.5±l GeV^/fm WM . 

Most of the empirical properties of the quenching factor for light-flavour hadrons - magnitude, 
Pj-, centrality-, ^/s^- dependences of the suppression - are in quantitative agreement with 
the predictions of non-Abelian parton energy loss models (Fig. [121). However, the fact that 
the high-pj- from semi-leptonic D and B decays is as suppressed as the light hadrons in 
central AuAu (Fig. [TTl right) [11391 [62l is in apparent conflict with the robust AEheavy-Q < 
^light-q < prediction of radiative energy loss models. Since the gluonsstrahlung 
probability of quarks and gluons is completely determined by the gauge structure (Casimir 
factors) of SU(3), the colour octet gluons (which fragment predominantly into light hadrons) 

II Technically, the q parameter can be identified with the coefficient in the exponential of an adjoint Wilson loop 



averaged over the medium length: (W^ (C) ) = exp (—1/4 ^/2)qL 



15J. 
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_ Suppression at pj = 4 GeV/c: 
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Figure 12. Left: Centrality dependence of Raa{pt > 7 GeV/c) for K° in AuAu and CuCu 
collisions at /s]^ = 200 GeV compared to the BDMPS-based PQM model lfT34l [1371 for 
different values of the (q) coefficient. Right: Excitation function of the nuclear modification 
factor, RAA{^/^), for %° in central reactions at a fixed = 4 GeV/c II138I . compared to 
predictions of a jet-quenching model with canonical non-AbeUan (soUd line) energy loss II141II . 



are expected to lose energy at Ca/Cf = 91 A times the rate of quarks. In addition, massive c, b 
quarks are expected to lose less energy than light ones due to their suppressed small-angle 
gluon radiation already in the vacuum ("dead-cone" effect) H 14011 . In order to reproduce the 
high-j9r open charm/bottom suppression, jet quenching models require either initial gluon 
rapidity densities (dN^/dy ~ 3000) [I133II inconsistent with the total hadron multiplicities, 
dN^ jdy ^ l.SdNch/dr\\^=o |!138] or with the dN^ jdy needed to describe the quenched light 
hadron spectra, or they need a smaller relative contribution of B relative to D mesons than 
theoretically expected in the measured decay electron pj range [11361 . This discrepanc>Q may 
point to an additional contribution from elastic (i.e. non-radiative) energy loss [|142[|143l[ for 
heavy-quarks [[133[[ which was considered negligible so far Will . The unique possibility at 
the LHC to fully reconstruct jets [|144| . to tag them with prompt y [I145[[ or Z [[146[[ and to 
carry out detailed studies in the c, h quark sector [[147[ [148[[ will be very valuable to clarify 
the response of strongly interacting matter to fast heavy-quarks, and will provide accurate 
information on the transport properties of QCD matter [fT5l[T6l . 

7. Propagation of collective perturbations in QCD matter: Distorted di-jet correlations 

Full jet reconstruction in AA collisions with standard jet algorithms [[149[[ is unpractical at 
RHIC energies due to low cross-sections for high-.E'T' jets and the overwhelming background 
of soft particles in the underlying event (only above ~30 GeV are jets above the background). 
Instead, jet-like correlations at RHIC are conventionally measured on a statistical basis by 
selecting high-pr trigger particles and measuring the azimuthal (A(|) = (|) — <^trig) and rapidity 

f Note, however, that the theoretical and experimental control of the pp ^ D,B +X reference (Fig. |5] left) is 
not as good as for the light hadron spectra 111 361 . 
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(At| = T] —V[trig) distributions of associated\iZ.drons ipT,assoc < PT,trig) relative to the trigger: 

1 d^N„ 
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' pair 
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ytrig dA(^dAr\ 

Combinatorial background contributions, corrections for finite pair acceptance, and the 
superimposed effects of global azimuthal modulations (elliptic flow) are taken into account 
with different techniques [|150lll511ll52[|153l . In pp or dAu collisions, a dijet signal appears 
as two distinct back-to-back Gaussian-like peaks around A(|) = (near-side) and A^ = n (away- 
side) {dAu panel in Fig.[T3l). At variance with this standard dijet topology in the QCD vacuum, 
early STAR results for semihard jets in central AuAu reactions HI SOU showed a complete 
disappearance of the opposite side peak for 3 < Ptmssoc < 4 < PTjhg < 6 GeV/c while 
the near-side correlation remained unchanged (Fig. [131 leftmost panel). Such a monojet-like 
topology confirmed a jet-quenching picture where a 2 — > 2 hard scattering takes place near the 
surface of the system with the trigger parton being unaffected and the away-side parton losing 
energy while traversing a medium opaque to coloured probes. For rising prjrig, the away-side 
parton is seen to increasingly "punchthrough" the medium II150L although the azimuthally- 
opposite correlation strength is still significantly reduced compared to dAu (Fig. [13]). 
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Figure 13. Angular correlations of high-p^ charged hadron pairs measured in: 0-5% central 
AuAu events in various PT.trig and pr.assoc (GeV/c) ranges (left), and in dAu and 0-5%-central 
AuAu for fixed pT,trig = 8-15 GeV/c and two pT,assoc ranges l,154J (right). 



The estimated energy loss, Eqs. ([8]) and ©, of the quenched partons is very large, up to 
AEto„ ^ 3 GeV/fm for a 10 GeV parton HI 371 and most of the (mini)jets, apart from those 
close to the surface, dump a significant fraction of their energy and momentum in a cell of 
about 1 fm^ in the rest frame of the medium. Since energy and momentum are conserved, the 
fragments of the quenched parton are either shifted to lower energy (pr <2 GeV/c) and/or 
scattered into a broadened angular distribution. Both softening and broadening of the away- 
side distribution are seen in the data [151 J when the px threshold of the away-side associated 
hadrons is lowered. Fig.[T4]shows the dihadron azimuthal correlations dNpair/dA^ in central 
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Figure 14. Azimuthal distributions of semihard hadrons (pr.assoc = 1 - 2.5 GeV/c) relative to a 
higher pr trigger hadron measured at RHIC. Left: PHENIX data in central (top) and peripheral 
(bottom) AuAu II 1521 (the arrows indicate the local maxima in the away-side hemisphere). 
Right: STAR data in central AuAu (squares) and dAu (circles) for two ranges of PT,trig II155I . 



AuAu collisions for PT,assoc = 1-2.5 GeV/c [|152[ I155 ll. In this semi-hard range, the away- 
side hemisphere shows a very unconventional angular distribution with a "dip" at A(|) ^ tt and 
two neighbouring local maxima at A(|) ^ n± 1.1. Such a non-Gaussian "volcano"-like profile 
has been explained as due to the preferential emission of energy from the quenched parton 
at a finite angle with respect to the jet axis. This could happen in a purely radiative energy 
loss scenario [I156II but more intriguing explanations for the conical-like pattern have been 
put forward based on the dissipation of the lost energy into a collective mode of the medium 
which generates a wake of lower energy gluons with Mach- [11571 11581 1159II or Cerenkov- 
like 11159111601 1161II angular emissions. In the Mach-cone scenario, the speed of soun^ of 
the traversed matter, = dP/de, can be determined from the characteristic supersonic angle 
of the emitted secondaries: cos(9m) = (c?), where 0^ is the Mach shock wave angle and (c,) 
the time-averaged value of the speed of sound of the medium traversed by the parton. The 
resulting preferential emission of secondary partons from the plasma at a fixed angle 9m ~ 1 • 1 , 
yields a value (cs) ~ 0.45 which is larger than that of a hadron-resonance gas (cs ~ 0.35) [I162II . 
and not far from that of a deconfined QGlEl (c^ = 1/ y^). Experimental confirmation of the 
Mach-cone picture or other alternative emission mechanisms for the associated particles in 
the quenched jet requires more detailed differential studies II163II such as the ongoing analysis 
of three- and many-particle azimuthal correlations II1641I165]| . 

+ The speed of sound is a simple proportionality constant relating the fluid pressure and energy density: P = e. 
* Note that although lattice calculations indicate that there are ^30% deviations from the ideal-gas limit in the 
s{T),P{T) and e{T) dependences up to very high T's [1|, the ideal-gas relation £ sa 3P (as well as other ratios of 
thermodynamical potentials) approximately holds above ~ 2T„i,. 
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8. Bulk hadronization: Enhanced baryon yields/flows at intermediate pT 

The increasingly suppressed production of mesons above pr ^ 2 GeV/c in central AuAu 
reactions at RHIC contrasts with the simultaneous unsuppressed p,p ni66l 11671 1168II and 
A, A i\69\ yields in the range pj "^l-A GeV/c. The intermediate-pr range mAuAu reactions 
features an anomalous baryon/meson~ 0.8 ratio which is roughly four times higher than in 
more elementary pp ox e+e^ interactions (Fig. [151 left)- Semihard (anti)protons show an 
enhancement with respect to the "xj- scaling" expectation for the ratio of perturbative cross 
sections at different cm. energies II170II . whereas xj scaling holds for all hadrons (mesons) 
measured mpp {AuAu) collisions at RHIC ifTFTlfTTOllTTTI . Only above ^ 6 GeV/c [[T68ll 
(Fig. [151 left) the meson/baryon ratio is again consistent with the expected yields obtained 
from universal fragmentation functions. Not only their spectra are enhanced, but at pr = 
2 GeV/c the V2 elliptic flow parameter of baryons exceeds that of mesons (v'2'^™" ^ 0.16, 
see Fig. |9l right) and keeps increasing up to pj- ~ 4 GeV/c when it finally saturates at 
^baryons ^ q 22 [|1721 I173II . All thosc observations clearly indicate that standard hadron 
production via (mini)jet fragmentation is not sufficient to explain the RHIC data for baryons 
at transverse momenta of a few GeV/c. 
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Figure 15. Left: Proton/pion in dAu and AuAu collisions at ^/s^ = 200 GeV 11681 compared 
to the ratio in light quark fragmentation in e+e^ at ^/s = 91.2 GeV (dotted-dashed line) and to 
parton coalescence predictions [41J . Right: Elliptic flow parameter V2 for all hadrons at RHIC 
normalised by the number of constituent quarks of each species (n^ = 2,3 for mesons,baryons) 
vs the transverse kinetic energy KEj = nij — m normalised also by 199111741 . 



The observed baryon nuclear modification factors close to unity, the high baryon/meson ratios, 
and their large elliptic flow have lent support to the existence of an extra mechanism of 
baryon production in AuAu based on quark coalescence in a dense partonic medium PTI . 
Recombination models compute the spectrum of hadrons as a convolution of Wigner 
functions with single parton thermal distributions leading to an exponential distribution 
which dominates over the standard power-law fragmentation regime below ~ 6 GeV/c. 
Coalescence models can thus reproduce the enhanced baryon production and predict also 
that the elliptic flow of any hadron species should follow the underlying partonic flow scaled 
by the number n of (recombined) constituent quarks in the hadron: V2{pt) = nv\{pT /n). 
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n = 2, 3 for mesons and baryons, respectively [11751 . Such a quark-scaling law for the elliptic 
flow is confirmed by the data [|172[ 11731 . Figure [TSl right, shows a recent variation of the 
quark-number scaling law for V2 that uses the transverse kinetic energy {KEt = mj — m, 
f^T = {Pt + m^y/^) rather than the pr and seems to account perfectly for the scaled V2 of 
all measured species also in the soft hydrodynamical regime below pr ~2 GeV/c fl99lll74ll . 
The overall success of valence quark coalescence models to explain hadron production in 
the semi-hard regime highlights the role of thermalised degrees of freedom in the produced 
system with partonic (as opposed to hadronic) quantum numbers. 

9. Temperature and equation-of-state (EoS): Thermal photons 

In order to describe the transient systems produced in AA collisions in terms of 
thermodynamical variables {T , £, s, etc.) linked by an EoS which can be compared to lattice 
QCD expectations, it is a prerequisite to establish that the underlying degrees of freedom 
form, at some stage of the reaction, a statistical ensemble. Proving that /ocaj^ thermalization 
has been attained in the course of the collision is thus a crucial issue both experimentally 
and theoretically [I176II . The large elliptic flow signal observed in the data strongly supports 
the idea of fast thermalization as discussed in Section [51 The identification of real and/or 
virtual y radiation from the produced "fireball" with properties consistent with a thermal 
distribution would in addition allow us to determine the underlying temperature and EoS of the 
system. Hydrodynamical [|40l and parton transport [il77l calculations indicate that the same 
cascade of secondary parton-parton collisions that drives the system towards equilibration in 
the first tenths of fm/c results in an identifiable emission of thermal radiation above the prompt 
perturbative yield in AA reactions in a window ^ 1 - 3 GeV/c. 
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Figure 16. Direct y,Y* spectra measured by PHENIX 11281 l96l in central AmAm collisions at 
y/sj^ = 200 GeV compared to NLO pQCD [64 1 and different hydrodynamics predictions |[92l 
including a QGP phase with Tq = 350-600 MeV. 

tt Note that, by simple causality arguments, global equilibrium in a finite system with radius Ra ^ 1 fm can only 
occur for time-scales x > 7 fm/c. 
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The preliminary y spectrum measured in central AuAu by PHENIX ll96l I178II is consistent 
with the sum of such perturbative (prompt) plus thermal (secondary) contributions (Fig.[T6l). 
Different hydrodynamical calculations of thermal photon production Il92lll79[|l80[|181[|182ll 
can reproduce the PHENIX data assuming the formation of a radiating QGP with tempera- 
tures around 2r„.„ {Tq = 400 - 600 MeV, corresponding to average temperatures (Tq) ^ 350 
MeV over the source profile) at times Xq = IRa/J^ 0.2 fm/c right after the two colliding nuclei 
pass-through each other. A word of caution should be noted in that the baseline (7^-scaled) 
prompt Y spectrum at = 1 - 4 GeV/c is obtained from NLO calculations, because no mea- 
surement is yet available in this pr range. The confirmation of the existence of a thermal 
enhancement over the prompt component, will require a direct measurement of the pp photon 
spectrum down to pj = ^ GeV/c. 

After subtracting the prompt y component from the inclusive direct photon spectrum, the local 
inverse slope parameter T^^ in the range pr ^ 2 - 4 GeV/c can be used as a relatively good 
surrogate of the initial medium temperature ^925. The combination of T.^ with another global 
observable directly related to the entropy density s would therefore allow one to determine 
the effective number of degrees of freedom g of the system via the Stefan-Boltzmann ratio 
g oc s/T^. Since the final charged particle rapidity density dNch/dr\\^=o is directly correlated 
with the initial entropy density of the system, by empirically studying the evolution of 
geff'^ dNch/dr\\r^=o/T^l versus T^g in different AA centralities one can effectively study the 
evolution of the number of degrees of freedom and look for any threshold behaviour related 
to the sudden increase at the transition temperature and/or a flattening of g.^ for temperatures 
above r„„. Such an approach has been tested in the context of a 2D-I-1 hydrodynamical model 
which effectively reproduces the hadron and photon data at RHIC ^2\. We found that one 
can clearly distinguish between the equation of state of a weakly interacting QGP and that of 
a system with hadron-resonance-gas-like EoS (i.e. with rapidly rising number of mass states 
with T). However, direct evidence of the parton-hadron phase change itself as a jump in g,ff 
around 7^^ ~ T„j, can only be potentially visible in AuAu reactions at lower center-of-mass 
energies ( ^/s^ ^ 20 - 65 GeV) II183H . At the LHC, the expected temperatures ^(1 GeV) 
reached in central PbPb will also produce a significant thermal photon signal up to pj ^ 6 
GeV/c [|184i which can be used to determine the thermodynamical conditions in the plateau 
regime of the s/T'^ EoS, closer to the ideal-gas limit than at RHIC. 

10. Critical temperature and energy density: Anomalous 7/\|/ suppression 

The study of heavy-quark bound states in high-energy AA collisions has been long since 
proposed as a sensitive probe of the thermodynamical properties of the produced medium [|37l . 
Analysis of quarkonia correlators and potentials in finite-T lattice QCD indicate that the 
different charmonium and bottomonium states dissociate at temperatures for which the colour 
(Debye) screening radius of the medium falls below their corresponding QQ binding radius. 
Recent lattice analyses of the quarkonia spectral functions L I 85 J indicate that the ground states 
(7/\|/ and T) survive at least up to T ^ 2r„,, whereas the less bounded %c and \|/' melt near 7^,„,. 
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Experimental confirmation of such a threshold-like dissociation pattern would provide a direct 
means to determine the transition temperature reached in the system and their comparison 
to ab initio lattice QCD predictions. A significant amount of experimental data on 7/\|/ 
production in different p{d)A andAA collisions has been collected at SPS HI 86[ 11871 II 8811 and 
RHIC [I189II . The corresponding nuclear modification factors compiled in HI 9011 are shown in 
Fig.[17]as a function of Npart- The surprisingly similar amount of 7/\|/ suppression observed 
at SPS and RHIC energies (with expected temperature differences of a factor of ~2) has been 
interpreted in a sequential-dissociation scenario [191J where the 7/\|/ survives up to T 2r„.„ 
in agreement with the lattice predictions, and the observed suppression at both cm. energies 
is just due to the absence of (30% and 10%) feed-down decay contributions from x{lP) and 
\|/'(25') resonances which melt already at T ^ r„„. The confirmation of such an interpretation 
would set an upper limit of T < 2r„„ ^ 400 MeV for the temperatures reached at RHIC. 
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Figure 17. nuclear modification factor versus centrality II 19011 (given by the number of 
participant nucleons in the collision) measured vciAA and p{d)A collisions at SPS 1, 1 861 1 1871 
[T88l and RHIC ll7T][T89l. 



Other explanations of the comparatively low depletion of 7/\j/ yields at RHIC have been 
put forward based on a much stronger direct 7/\|/ suppression (at temperatures close to 
7^.,,,) combined with cc pairs regeneration from the abundant charm quarkqijj in the dense 
medium II192II . The LHC measurements will be crucial to resolve this issue. A strongly 
suppressed 7/\[/ yield in PbPb at 5.5 TeV - where the expected initial temperatures will be 
well above 27^,„, - would support the sequential- screening scenario, whereas recombination 
models predict a strong enhancement due to the larger density of cc pairs in the medium. In 
addition, the abundant production of the T{ls,2s,3s) states at LHC energies will open up a 
unique opportunity to study the threshold dissociation behaviour of the whole bottomonium 
family. The T is expected to survive up to 47^,,,, and, therefore, direct suppression of the bb 
ground- state would be indicative of medium temperatures around 1 GeV at the LHC. 



ttlO charm pairs are produced on average in a central AuAu collision at the top RHIC energy. 
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In-medium modifications of the spectral function (mass, width) of the light vector mesons 
(p, (0, and (|)) have been proposed as a promising signature of the (approximate) restoration 
of chiral symmetry in the u,d,s quark sector [|38l [391 . In the QCD vacuum the spontaneous 
breaking of chiral symmetry manifests itself in the hadron spectrum through the mass- splitting 
of "chiral partners" i.e. between states of opposite parity but equal quantum numbers. Chiral 
symmetry breaking leads, in the mesonic sector, to the non-degeneracy of the pseudo-Zscalar 
(tt — a) and axial-Zvector {a\ — p) channels. For temperatures above the chiral transition, 
massless left- and right-handed quarks will decouple and one expects to observe a gradual 
disappearance of the mass-splitting, leading to a shift of the masses of vector mesons and 
their chiral partners. The p meson is an excellent candidate for the experimental study of in- 
medium spectral functions in AA collisions due to (i) its short lifetime (Xq = 1 .3 fm/c) allowing 
it to decay before regaining its vacuum spectral shape, and (ii) its (rare but detectable) dilepton 
decay branching ratio (r/_|_/_ ^ 5 ■ 10~^) which is unaffected by final-state interactions with 
the surrounding environment. 
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Figure 18. "Excess" invariant /j^/j^ mass spectrum in semi-central Inin collisions at 
y/s^= 17.3 GeV [32] compared to the expected p line-shape in different theoretical scenarios: 
vacuum p, collisional broadening Il38lll93l . and mass drop ||39l . 

The NA60 experiment has recently studied a high-statistics sample of low mass muon pairs 
in Inin collisions at y/s^= 17.3 GeV [|32l . The measured p spectrum, after subtraction of all 
other sources of opposite-sign muon pairs, is peaked at the nominal (free) mass (mp = 0.77 
GeV/c^) but its width is substantially broadened for increasing centralities (Fig. [TS]) . The p 
spectral shape is better reproduced by models which consider in-medium width broadening 
due to interactions in a hot and dense hadronic medium close to the expected phase boundary 
(r ~ 190 MeV) II193II , than by a long-standing prediction based on a downwards mass- shift 
coupled directly to the melting of the chiral condensate [|39ll . Such a result sets new constraints 
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on the possible realization of chiral symmetry in QCD matter. At RHIC energies, the factor 
of two larger temperatures reached compared to SPS prefigure that the p should exhibit a 
completely "melted" line-shape if, as expected from the lattice, the chiral and deconfinement 
transitions occur at the same 7^,;,. The p measurement in the dielectron channel will be 
possible in PHENIX with the recently installed hadron-blind-detector HI 9411 which will allow 
to suppress by two orders of magnitude the large combinatorial background arising from light- 
meson Dalitz decays and photon conversions, while preserving 50% of the signal. 

Summary 

High-energy collisions of heavy ions provide the only existing method today to explore em- 
pirically the phase diagram of QCD at extreme values of temperature, density and low-jc. 
We have reviewed the experimental and phenomenological progress in nucleus-nucleus col- 
lisions at BNL-RHIC ( ^ 20 - 200 GeV) and CERN-SPS ( ~ 20 GeV) in recent 
years with a particular emphasis on the modifications suffered by different hard probes in 
the QCD medium that is produced, compared to baseline measurements in proton-proton or 
proton(deuteron)-nucleus collisions. The observed modifications allow us to obtain direct in- 
formation on fundamental (thermo)dynamical properties of strongly interacting matter. 

The reduced AuAu hadron multiplicities and the depleted yields of semihard hadrons {pT ~ 
1-4 GeV/c) at forward rapidities in dAu collisions indicate that the initial conditions of the 
nuclei accelerated at RHIC energies are consistent with Colour-Glass-Condensate approaches 
that model the hadronic parton distributions as a low-x saturated gluon wavefunction evolving 
according to non-linear QCD equations. These new data - complementary to the rich HERA 
results on the partonic structure of the proton - shed new light on the high-energy limit of 
QCD, a physics topic not only appealing in its own right but an essential ingredient for any 
serious attempt to compute a large variety of hadron-, photon- and neutrino- scattering cross 
sections at increasingly large energies. 

The robust radial and elliptic flows seen for all identified hadron species up to /Jy- ~ 2 GeV/c 
in AuAu reactions at ^/s^ = 200 GeV are remarkably well described by ideal relativistic hy- 
drodynamics calculations which model the expanding system starting with a realistic QGP 
equation-of- state with initial energy densities Eq ~ 30 GeV/fm^ at thermalization times Xq ^ 
0.6 fm/c. Detailed hydro-data comparisons for various differential observables - in particular 
the elliptic flow parameter V2{pt) for different light and heavy hadron species - indicate that 
the system exhibits very small viscosities (i.e. very short mean free paths) and is strongly cou- 
pled at variance with the anticipated QGP paradigm of a weakly interacting gas of relativistic 
partons. In the intermediate pr ^2-4 GeV/c range, baryons have enhanced yields and flows 
compared to mesons pointing to a novel channel for hadronization based on constituent-quark 
coalescence in a dense partonic medium. The existence of a new mechanism of hadron forma- 
tion in heavy-ion reactions at transverse momenta of a few GeV/c, apart from standard parton 
fragmentation, offers new ways to probe the space-time dynamics of confinement in different 
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QCD environments. 

In the high-/?r sector, leading hadrons (but not colour-blind prompt y) are suppressed by 
up to a factor of ~5 in the range f^i 4 - 20 GeV/c compared to perturbatively-scaled 
proton-proton spectra. This result is in agreement with non-Abelian energy loss models that 
assume that the produced partons traverse a medium with very large parton rapidity densi- 
ties dN^/dy 1100 and transport coefficients (4) ~ 10 GeV^/fm. The energy lost by the 
quenched parton in the medium apparently shows up in a very unconventional conical-like 
azimuthal profile of secondary hadrons (pj ~ 1 - 3 GeV/c) in the away-side hemisphere of 
high-pT trigger hadrons. Interpretations of this preferential azimuthal emission at A(|) 7t± 
1.1, as caused by the generation of a Mach-cone boom by a supersonic parton propagating 
through the dense system, yield average speeds of sound (c^) = cos (0m) ~ 0.45 not far from 
those expected from lattice QCD for deconfined quark-gluon matter. 

In the electromagnetic sector, preliminary real and virtual y spectra in central AuAu in the 
range = 1 - 14 GeV/c can be described as the sum of a perturbative (prompt) photon 
contribution plus a secondary component of thermal origin. Hydrodynamics calculations can 
reproduce the data assuming the formation of a radiating QGP with average temperatures 
2 r^„, 350 MeV. Such temperature values seem to be consistent with the observation that the 
amount of J/y\f suppression at RHIC and SPS is about the same as expected from lattice-based 
calculations that predict a survival of the cc ground state up to T ^ 2r„„ but a "melting" of the 
Xc and V/' (which feed-down at a ~40% level to the 7/\|/) near T„,,. At the CERN SPS, recent 
NA60 results on the p spectral function in central Inin collisions at ^/s^ =17.3 GeV indicate 
that the width of the vector meson is substantially broadened in the medium. Theoretical cal- 
culations indicate that most of the broadening can be accounted for by coUisional effects in a 
hot and dense hadronic medium with initial temperatures close to T^^, ftil75 MeV (note that 
such a result is also consistent with the observed hadron abundances at SPS which indicate 
that the system reaches chemical equilibrium at temperatures around 160 MeV). 

The overall scenario taking form from the wealth of recent experimental data suggests, on the 
one hand, that heavy-ion collisions at RHIC energies produce a strongly interacting liquid-like 
QGP with very large initial parton rapidity densities dN^/dy ~ 1100, temperatures 2r„„ ^ 
350 MeV and very low shear viscosities. On the other hand, systems produced at SPS seem 
to be only partially thermalised - according to the lower measured collective anisotropic flow 
compared to RHIC -, have initial dN^/dy ^ 400 and temperatures around the phase boundary 
at T^rtt ~ 175 MeV. However, it is fair to acknowledge that among the existing signals there 
is yet no incontrovertible "textbook" figure proving the formation of a thermalised extended 
medium consisting of deconfined and chirally-symmetric quarks and gluons. One such evi- 
dence would be a direct empirical observation of a jump in the number of effective degrees 
of freedom at the phase change as expected by EoS calculations in the lattice. Since RHIC 
top energies seem to produce systems at twice r„„ and SPS data point to conditions just at 
the predicted phase change, the expected jump is likely to be observed (e.g. via precise stud- 
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ies of the correlation of thermal photon slopes with the global hadron multiplicities) in a next 
phase of RHIC running at intermediate energies ( ^/s^ ^20-62 GeV) and high luminosities. 
Likewise, confirmation of the concurrent chiral and deconfinement transitions in QCD matter 
will require e.g. precise measurements of the p spectral shape (and, ideally, that of its chiral 
partner a\) at RHIC energies. Lower-energy runs at RHIC, as well as at the projected CBM 
facility [| 19511 , will access also the region of large baryon densities around the QCD critical 
point. Finding signs of the tri-critical point at relatively high temperatures would indicate that 
the smooth cross-over changes to a first order phase transition at higher baryon densities, a re- 
sult of relevance for the conditions prevailing in the core of neutron (and other compact) stars. 
Direct validation of the strongly-coupled interpretation of the medium formed at RHIC and 
potential observation of the anticipated weakly interacting quark-gluon plasma require key 
measurements in PbPb at 5.5 TeV at CERN-LHC where the initial temperatures ^(1 GeV) 
should be large enough to observe the direct melting of ground-state quarkonia resonances. 
In addition, at the LHC, the longer duration of the QGP phase and the much abundant pro- 
duction of other hard probes (in particular parton energy loss results for fully reconstructed, 
y- or Z-tagged, and flavour-identified jets) thermal photons, V2 flow parameter, etc. will likely 
result in indisputable probes of the deconfined medium much less dependent on details of the 
later hadronic phase. 

The experimental advances in the last years have been paralleled by significant progresses 
in the theoretical description of high-density QCD matter. Lattice methods are increasingly 
more refined and powerful to describe not only the infrared collective dynamics (EoS, criti- 
cal parameters) but also the in-medium correlators (quarkonia). Effective field theories have 
been developed in specific domains such as the Colour-Glass-Condensate which effectively 
describes the high-energy (low-x) limit of QCD. Perturbative calculations have substantially 
improved the description of the interaction of hard probes with hot and dense quark-gluon 
matter. Last but not least, duality approaches based on the application of AdS/CFT cor- 
respondence between weakly coupled gravity and strongly coupled QCD-like systems, are 
providing new powerful insights on dynamical properties that cannot be directly treated by 
either perturbation theory or lattice methods while simultaneously opening novel directions 
for phenomenological studies and experimental searches. 

The impressive experimental and theoretical advances triggered by the wealth of high- 
statistics, high-quality data collected in ultrarelativistic nucleus-nucleus collisions at RHIC 
and SPS, have significantly expanded the knowledge of many -body QCD at extreme 
conditions of temperature, density and low-.x;. Those studies - which will be substantially 
extended in the upcoming LHC (and likely RHIC-II) nucleus-nucleus and proton-nucleus 
programme - go beyond the strict realm of the strong interaction and shed light on a vast 
ramification of fundamental physics problems. Knowledge of the collective behaviour of 
many-parton systems is of primary importance not only to address basic aspects of the strong 
interaction such as the nature of confinement or the mechanism of mass generation via chiral 
symmetry breaking, but to ascertain the high-energy limit of all scattering cross sections 



Quark-Gluon Matter 



28 



involving hadronic objects, the inner structure of compact stellar objects, or the evolution 
of the early universe between the electroweak transition and primordial nucleosynthesis. 
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